Nonlinear graphene plasmonics: amplitude equation 
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Using perturbation expansion of Maxwell equations, the amplitude equation is derived for non- 
linear TM and TE surface plasmon waves supported by graphene. The equation describes inter- 
play between in-plane beam diffraction and nonlinerity due to light intensity induced corrections 
to graphene conductivity and susceptibility of dielectrics. For strongly localized TM plasmons, 
graphene is found to bring the superior contribution to the overall nonlinearity. In contrast, non- 
linear response of the substrate and cladding dielectrics can become dominant for weakly localized 
TE plasmons. 
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I. INTRODUCTION 

Applications of graphene in photonics and optoelec- 
tronics are being actively discussed in recent years P, 0] ■ 
In particular, graphene plasmonics is considered as a 
promissing alternative to conventional plasmonics with 
noble metals Q. Recently, hybrid metal-graphene plas- 
monic structures have been proposed as the propitious 
platform for novel optical devices Q. 

Graphene supports two types of surface plasmons: 
transverse magnetic (TM) and transverse electric (TE) 
modes @|. TM graphene plasmon is in many 
ways analogous to the surface plasmon excited at a 
metal/dielectric interface Q, although specific features of 
collective electron excitation in the purely 2D graphene 
lead to qualitative differences in the spectra of plasmons 
in these two systems @. Compared to its metal analogue, 
TM plasmon supported by graphene offers susbstantial 
enhancement of the field localization, accompained by 
the considerable decrease of the propagation loss - all be- 
ing crucial for potential applications of surface plasmons 
in miniature photonic components. The existence of TE 
plasmon is directly related to the linear (Dirac) spectrum 
of electrons in graphene |5| , there is no analogue of such 
surface wave in conventional plasmonics. TE plasmon is 
only weakly localized at the surface, however it is char- 
acterized by considerably low propagation losses even at 
room temperatures. Spectral characteristics of TE and 
TM graphene plasmons are defined by the charge density, 
which can be controlled chemically Q or electrically (9j. 
This tunability represents another important advantage 
of graphene plasmons over metal plasmons. 

Optical properties of doped graphene are encapsulated 
in the induced surface current /C. So far, graphene plas- 
mons have been studied under the assumption of a linear 
relation between the current K, and the field amplitude 
E: JC = o£ 19]. This is true only at low light intensities, 
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FIG. 1: (Color online) Schematic illustration of surface plas- 
mon propagating along graphene sheet. Fields are exponen- 
tionally localized in x (across the interface), as shown in the 
left panel. 



servation of four-wave mixing with graphene flakes [14j 
and a range of nonlinear effects in a graphene-coated pho- 



tonic crystal nano-cavity 15j. Altogether these findings 



while generally the dependence K{£) is predicted to be 
highly nonlinear [Tol \\A\ . The particularly strong nonlin- 
ear repsonse of graphene has been confirmed in several 
experiments, including direct measurements with optical 
Kerr gate [l2| and z-scan [l3| techniques, as well as ob- 



put forward the great potential of graphene for building 
functional nano-photonic devices. 

In this work we consider nonlinear surface waves sup- 
ported by graphene in the simple planar geometry shown 
in Fig. [TJ Allowing for light intensity corrections to the 
surface current and to the susceptibility of dielectrics sur- 
rounding graphene, as well as introducing diffraction due 
to a finite beam width in the unbound (y) direction, we 
develop asymptotic expansion of Maxwell equations and 
boundary conditions to obtain an amplitude equation for 
quasi-TM and quasi-TE surface waves. The asymptotic 
expansion procedure is similar to that recently developed 
for semiconductor and metal nano- waveguides [H, u!f[ . 
Further, we analyze the relative contribution from di- 
electrics and graphene to the overall effective nonlinear- 
ity of the system for the two types of plasmons, and the 
impact of geometry on the nonlinearity enhancement. 



2 



II. SETUP AND ASYMPTOTIC EXPANSION 
OF MAXWELL EQUATIONS 

Wc consider the planar geometry, in which single layer 
graphene is sandwiched in-between two dielectrics. We 
choose x axis to be perpendicular to the interfaces, z is 
the direction of propagation, and y is the unbound direc- 
tion in which light can diffract, see Fig.[TJ For monochro- 
matic fields, £ = \Ee-' lbJt + c.c, H = \He- iult + c.c, in 
each dielectric domain we solve stationary Maxwell equa- 
tions: 



V x V x E = — 



(1) 



Here spatial coordinates are normalized to the inverse 
wave number k — 2tt/X = uj/c. For homogeneous 
isotropic dielectrics, dispacement vector takes the form: 



D = e eE 
1 



N 



:X3 



+ N 
E\ 2 E 



-E 2 E* 



(2) 
(3) 



Adapting complex amplitude notation to the surface 

the 



current K. = \Ke~ lult 



c.c, K = [0,K y ,K z 
boundary conditions can be written as: 



A[E y ] = , A[E Z ] = , 
-A[H V ] = ice A[d z E x - d x E z ] = K z , 
-A[H Z ] = ice a A[d x E v - d y E x ] = -K y , 

where operators A and are defined as: 

A[f{x)\ = lim(f(x-S)-f(x + 6)) , 



©[/(*)] 



lim (f(x-S) + f(x + S)) 

2 5->0 



(4) 
(5) 
(G) 



(7) 
(8) 



Taking into accound first order nonlinear corrections to 
the relationship betwen surface current and electric field 
amplitude JC{8) [ll[ and neglecting the effect of higher 
harmonics generation, we obtain: 



highly doped graphene, ^> kT and Tllo < 2\/j,\, /j, is the 
chemical potential [5]. Furthermore, we assume that non- 
linear corrections to the dielectric susceptibility ~ X3l^| 2 
and graphene conductivity ~ cr 3 |i?| 2 are of the same or- 
der of smallness O(s). We let the mode amplitude tp to 
vary slowly with propagation distane, d z tp <C /3ip, and 
consider weak diffraction, d y ip ^ 0. Using asymptotic 
expansion of Maxwell Eqs. (JT|) and boundary conditions, 
Eqs. ((I])-©, below we derive propagation equation for 
the mode amplitude tp. 

Note, with the account of diffraction, the separation 
into TM and TE modes can no longer be performed, in- 
stead one deals with quasi-TM and quasi-TE modes. 



III. QUASI-TM SURFACE PLASMON 

We seek a guided mode solution in the form: 



E , 



A x (iP,x) + B x (ip,x) + 0(s 5 / 2 ) 



E y = [C(tP,x)+0(s 2 )]e i ^ , 



(10) 
(11) 



A z {i>,x)+B z (iP,x) + 0(s 5 / 2 ) e^ 2 , (12) 



where ip — ip(z, y) is a slowly varying function: d z tp ~ s, 
d y ip ~ s 1 ' 2 , A x>z ~ s 1 ' 2 , C ~ s, B ~ s 3 / 2 . The chosen 
orders of smallness are justified below by solving con- 
sistently boundary value problems, that emerge in dif- 
ferent orders of s. Following substitution into Maxwell 
equations, in the order O^s 1 ^ 2 ) we obtain the following 
boundary value problem: 

LtmA = , (13) 
A[A Z ] = , A[ipA x - d x A z ] = o^ ] Q[A z ] , (14) 

where A = [A x , A Z ] T , a% = a\/{ceo), and operator Ltm 
is defined as: 



Ltm = 



ft 2 



i/3d x 



i/3d x -d 2 



(15) 



K y , z = e 
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2 



E\ 2 E 



1 



y,z + 7> E & 



y.z 



(9) 



If the nonlinear response is neglected altogether (er 3 = 
0, X3 = 0)i an d n0 losses due to electron-phonon scat- 
tering or defects are considered at zero temperature 
(Re(ai) = 0), the above system admits solutions in 
the form of surface plasmons propagating in z direc- 
tion: E,H ~ e l P z . Field amplitudes in such solutions 
are exponentially localized at the interface x = and 
constant along the unbound direction y. For the case of 
positive/negative imaginary part of conductivity o\ only 
TM/TE surface plasmon exists [5[. 

Below we consider a weakly dissipative case, so that 
a i = a[ R ' + ierj 7 ' and a[ R ^ /a[^ ~s<l, where s is a 
dummy small parameter. This assumption is valid for a 



We choose the solution in the form A — I 1 / 2 ip{z,y)e, 
where e = [e x , e z ] T is the linear surface plasmon mode: 



x > 



-ip 



i(3 



(16) 

(17) 
(18) 



e s and e c correspond to dielectric layers at x < (sub- 
strate) and x > (cladding), respectively. Propagation 
constant /3 is defined through the dispersion relation Q : 



(19) 
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The normalization factor / is chosen in a way that 
\ip\ 2 is the power density (measured in watts per meter) 
carried in the z direction [la ]: 



/ 

Q 



2(3k 
e cQ 

+ 0O 



eleJ dx 



/3 2 



(20) 
(21) 



2 <K 
Collecting terms of the order O(s) we obtain: 

((3 2 - e)C - d 2 xx C = -I 1/2 d v i; (ipe z + d x e x ) , (22) 
A[C] = , A[d x C - I^ 2 d y ipe x ] = -a{ I} e[C} .(23) 

From V • E = in the order 0(s 1/2 ) it follows that 
if3e z + d x e x = 0, and therefore C solves the homogeneous 
equation. It is non-zero due to simultaneous diffraction 
[dyi]) ^ 0) and discontinuity of e x component at the in- 
terface, see Eq. (|23|) . Substituting C — E^dyipey, it is 
easy to see that e y satisfies the same homogeneous equa- 
tion as e z . Comparing boundary conditions for e y and 
e z , we obtain e y = (—i//3)e z . 

In the order 0(s 3 / 2 ) we obtain the following boundary 
value problem: 



LtmB — I 1 / 2 J , 
A[B Z ] = , 



(24) 
(25) 



A[ipB x + d z ipi 1 / 2 e x -d x B z ] 



+a[ J) e[B z ] - ia 3 / 3/2 |^|Ve [\e\ 2 e z 
where a 3 — (73/(060) and 

J x = d z i>{2i/3e x - d x e z ) + d 2 y ip(e x - d x e y ) 

+|V'| 2 V' n x , 
Jz = -d z ipd x e x + d 2 y ip(e z - i/3e y ) 

+ |V>| 2 V' n z , 



(26) 



(27) 
(28) 



n x , z = IN x , z (e x ,0,e z ). 

Next, we project Eq. (IM|) onto the linear mode e: 



(e* • L TM B} dx = I 1/2 f + (V ■ fj dx , (29) 



take J_ = J_ + L in the l.h.s., apply integration 
by parts and use boundary conditions in Eqs. (1141) and 
(1561) to obtain: 



dx 



J + (e* • UmB 

A [i/3(e* z B x + e* x B z ) - e* z d x B z + B z 8 x e* z ] 

= -d z ^ 2 A[e x e* z ] - ia[ R) 4>I 1,2 Q[\e z \ 2 ] 

1 



^a 3 |^|V/ 3/2 e 



|eT|e z | 2 + ^(e:) 2 



(30) 



Finally, computing integrals in the r.h.s. of Eq. f|29[) . 
we obtain the amplitude equation: 

<hr +?AV; + T|^| 2 ^ = 0, (31) 



dip 



+ 



1 



d(z/k) 2/3kd(y/k) 2 

where the nonlinear parameter T combines contributions 
of graphene and dielectrics: 



T 

7G 

P 



g{iG + iD) 

■ia 3 k 2 
2e c/3 2 P 2 " 
k 2 j 
2e cf3 2 P 2 J_ 



|ei 2 | ez | 2 + ^ 2 ( e ;) 2 



+00 



+00 



\e\ 2 dx = 



X3 
2/3 2 - e t 

2gf 



2 16 



2/? 2 



2|2 



2gl 



(32) 
, (33) 

, (34) 

■ , (35) 



the surface-induced nonlinearity enhancement factor g is 



(! + »?)■ 



T 1 + 1 

and the effective linear absorption parameter is given by: 



(36) 
(37) 



.,(«) 



A 



9 



1/2"! 



2/3P 



e[|e,| 2 ] 



(38) 



In the above derivations we used the auxiliraly relation 
g x / 2 Q = j3 2 P, which can be obtained by usin g %l3 e r . = 
—d x e x and taking by parts integral in Eq. (|55|) [16]. 

Expression for graphene nonlinear coefficient in Eq. 
(|33f can be replaced by the integral similar to the one in 
Eq. ([M]), following introduction of an effective graphene 
nonlinear susceptibility: 



xi r) = ia 3 S(x) 



r (3) 



e c 



-S(x) 



(39) 



where S(x) is the Dirac delta function. Note however 
the different structure of the term under the integral, 
which is due to the surface nature of nonlinear response 
in graphene, cf. Eqs. (JSJ) and In the limit of high 
localization, f3 3> e SlC , for the guided mode one obtains 
simple relation e x — He z , and therefore: 







(40) 



Apparently, in this limit, the effective nonlinear response 
of graphene is twice weaker than that of a infinitesimally 
thin Kerr medium with the susceptibility xi 



(a r ) 



IV. QUASI-TE SURFACE PLASMON 

For the case of quasi- TE mode we use the ansatz: 
E x = [C x (il>,x) + 0(s 2 )]e^ z , (41) 



Ey - 

E, = 



Adj } x)+B^j,x)+0( S 5/2 ) 
[C z ^,x)+0(s 2 )] e*' , 



J0Z 



(42) 
(43) 
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where d z ip ~ s, dyip ~ s 1 ^ 2 , A ~ s 1 / 2 , C x . 2 ~ s, S ~ 
s 3 / 2 . Following substitution into Maxwell's equations, in 
the order 0(s 1 ^ 2 ) we obtain the following boundary value 
problem: 



L TE A = 0, 
A[A] = 0, A[d x A] = -a[ I) e[A] 
Lte = P 2 - 6 ~ dlx ■ 



(44) 
(45) 
(46) 



We choose the solution in the form A = J 1 / 2 -0(2, y)e 
where e y is the surface plasmon mode: 

X < : e y = e qsX , 
x > : e„ = e~ 9ca: , 



(47) 
(48) 



q s . c arc defined in Eq. (|18p. and the normalization factor 
/ ensures |^| 2 gives the power density carried in the z 
direction: 



P 



4A- 



+ 00 



'da; 



1 
2q s 



1 

2fc 



Dispersion relation for the TE plasmon is given by: 



VP 



VP 2 



-a 



(I) 



(49) 
(50) 

(51) 
(52) 

A[C Z ] = , A[0C x - d x C z ] = a^ ] Q[C z ] , (53) 

where C — [C X ,C Z ] T , operator Ltm is defined in 
Eq. (fT5|). Substituting in the above equations C — 
l x l 2 d y ^}e and eliminating e x , we obtain: 



In the order O(s) we obtain: 

LtmC = -I 1/2 dyij[d x e y ,if3e y ] T 



LteZz = ifie 1 L TE e y = . 



(54) 



In other words, e z solves the same homogeneous equation 
as e y . Comparing boundary conditions for e y , Eq. ([4"5"j). 
and e z , Eq. (|53|) . we obtain e z = (i/ f3)e y and e x = 0. It is 
easy to check that this choice also satisfies the condition 
V • E = in the order O(s). 

In the order 0(s 3 / 2 ) the following boundary value 
problem is obtained: 

L TE B = I 1 ' 2 [i2pd z ^e y - d 2 y iP(i(3e z + 3 x e x )] ,(55) 
A[B] = , (56) 
A[d x B - d y C x ] = ia^lV^Q^] 

-a^SlB] + t la 3 P/ 2 \^\ 2 m\e y \ 2 e y } ,(57) 

Projecting Eq. ([SB"]) onto the mode e y and following essen- 
tially the same steps as described in the previous section, 
in the l.h.s. we obtain: 



«+oo 



e*LTEBdx 



-Q[e*d x B - Bd x e*} 



(58) 



Performing projection in the r.h.s. of Eq. (|55[) . and using 
boundary conditions in Eqs. (|45j) and (|57j) . we obtain the 
amplitude equation ([3Tj) with the following coefficients: 



1G,TE 
1D,TE 



1G,TE + 1D,TE ■ 

i3a 3 k 2 r 
2e cP 2 P 2 [l6y 
3fc 2 



A 



TE 



2e c(3 2 P 2 



a [R) k 

^e\\e y \ 2 } 



X3\e y \ 4 dx , 



2/3P 



(59) 
(60) 

(61) 
(62) 



For quasi-TE mode the enhancement factor g is absent. 
Also, due to the linear mode being scalar, in this case the 
nonlinear response of graphene is completely analogoues 
to that of a infinitesimally thin Kerr medium with the 
susceptibility x 3 sr) in Ec l- <E3>- 



V. ANALYSIS AND DISCUSSION 

Conductivity of graphene consists of intra- and inter- 
band contributions, o\ = <J intra + cr inter . For the case 
of a highly doped graphene, \fj,\ 3> kT, /i is the chemical 
potential, intra- and inter-band terms are given by the 
semi-classical formalism Q : 



&intra (^) 
Winter (^) 



ie 2 1 

irk Q + iVintra 
ie 2 2 — \fl\- iv in ter 

In j — j 

Anh 2+ O \+lVinter 



(63) 
(64) 



where Q = hw/fjt, excitation below interband absorption 
threshould is assumed: 57 < 2, coefficients v = ?i/(|/x|r) 
take into account losses due to electron scatterings at 
finite temperatures, below we take Ti ntra = lOOfs and 
Unter = lps @, HH ■ For the doping level of fi = O.leV 
we obtain v intra 0.066, v inter w 0.007, and the inter- 
band absorption threshould is at U)th = 2/i/h s» 3 • 10 14 
rad/s (Xth ~ 6.3/im). The corresponding dimensionless 
conductivity a\ is plotted in Fig. [21 Imaginary part of 
cti changes its sign at Slo ~ 1.67, linear TM plasmons 
exist for £1 < fig (i.e. when > 0), while TE - for 

> Qo («i < 0) 5]. For the chosen doping level, f2 
corresponds to Ao ~ 7.5/xm. 

The nonlinear conductivity coefficient 03 for graphene 
is given by [Til ]: 



2t2 



0-3 (fi) = 



(eV F ) 2 h 



32 irh /i 4 r2 3 



-(1 + iar) ■ 



(65) 



Here we introduced coefficient ar to account for two pho- 
ton absorption in graphene, recent experiments suggest 
a T sa 0.1 I5j- Negative imaginary part of er 3 suggests 
that the nonlinear response is of self-focusing type, cf. 
Eqs. (EJ), (021) and (UHl). 
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FIG. 2: (Color online) Dimensionless conductivity of 
graphene ai: imaginary (a) and real (b) parts. Dashed, dot- 
ted and solid curves correspond to intra-, inter- and full con- 
ductivity, respectively. The chemical potential is set to /i — 
O.leV, and relaxation times are Ti„tra ~ lOOfs, Tinter ~ lps. 

Below we consider surface plasmons in configurations 
with air (e = 1) and silicon (e = 12) as dielectrics. For 
silicon we take %3 = (4/3)ceoe s ri2, ni = 4 • 10~ 18 m 2 /W. 
Two photon absorption in silicon is negligible for A > 
2/im |18( . For simplicity we neglect dispersion of linear 
and nonlinear dielectric constants. 



A. Quasi-TM plasmon 

First, we consider TM plasmon in the configuration 
with silicon substrate and air cladding. The correspond- 
ing dispersion is plotted in Fig. [3] To validate our the- 
ory, the propagation loss parameter A/k is compared 
against imaginary part of the propagation constant (see 
open squares in Fig. [3]) . The latter is computed from the 
full dispersion relation that takes into account complex- 
valued «i and obtained by replacing with — ia\ in 
Eq. (jTHJ) . The results are found to be in perfect agree- 
ment. 

TM plasmon is characterized by a considerable light 
confinement in a wide range of frequencies: even at the 
frequency as low as ft — 0.25 (lu — 7.5 • 10 13 rad/s, A ~ 
25/im) the propagation constant is (3 w 240, and it con- 



OQ. 




FIG. 3: (Color online) Dispersion of TM plasmon in the 
silicon-graphene-air geometry: propagation constant f3 and 
loss parameter A/k as functions of Q. Imaginary part of f3 
as computed from the full dispersion relation with complex- 
valued qi is shown with squares. 



stantly grows as increases towards the threshould value 
£Iq. Propagation losses are relatively low: A/(kf3) < 10 -3 
when 0.5 < < 1.6, which is due to the low absorp- 
tion rate in graphene below the interband absorption 
threshould, cf. Fig.^b). 

Due to the high localization of TM plasmon, the 
surface-induced enhancement factor g is large and is 
growing nearly exponentially with increasing Q, see 
Fig. BJa). This growth over-ballances the decay of 
graphene nonlinear response 03 ~ £^~ 3 , cf. Eq. (|65l) . 
and causes the considerable increase of the nonlinear co- 
effficient T with frequency, see Fig.HJb). 

Remarkably, the relative contribution of silicon sub- 
strate to the overal nonlinearity remains negligibly small 
within the entire frequency window of existence of TM 
plasmon, as illustrated in Fig. EJb). Also, due to the 
large j3, the diffraction term in Eq. (|3"TT) can be neglected 
for typical beam widths L y > Ifim. Indeed, for f2 = 1 
(A 12.4y«77i) and the beam width of L y = 10/im the 
diffraction length is Lq = L y (3k m 50mm, which is more 
than six orders of magnitude larger that the apparent 
plasmon wavelength X p = 2%/(pk) ~ lOnm. Neglecting 
nonlinear response of dielectrics and beam diffraction, 
as well as disregarding linear and nonlinear absorption 
(a[ R ^ — = 0), one can find stationary solutions 

of Maxwell equations E(x, y, z) = Iip e(x, y; (3NL)e l ^ NLZ 
with the nonlinear boundary condition in Eqs. (|4][6]) an- 
alytically. The corresponding dispersion relation reads: 

+ ^l[\^\e z \ 2 + \^{elf]\M 2 -(66) 

At the same time, substituting i/j(y,z) = ij]Qe l ^ NLZ into 
the amplitude equation (|31[) and assuming A = Im(T) = 
0, we obtain Pnl = ft + {T / k)\^ \ 2 . Solving the dis- 
persion relation in Eq. (1661) numerically, we found both 
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FIG. 4: (Color online) Effective nonlinearity of quasi- TM 
plasmon: (a) surface enhancement factor g; (b) nonlinear co- 
efficient T and the relative dielectric nonlinearity. 




FIG. 6: (Color online) Dispersion of TE plasmon in the 
silicon- gr aphene- silicon geometry: Plasmon localization fac- 
tor q — \/ P 2 — e, propagation loss A/k and imaginary part of 
the propagation constant (squares). 
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FIG. 5: (Color online) Nonlinear index shift A/3 = /3nl — ft 
vs power density. Solid curves correspond to the numerical 
solution of the dispersion relation in Eq. (I66I) . dashed lines - 



to the result given by the amplitude Eq. ()31[) : A/3k = T | "J^q 



results to be in perfect agreement at low power densities 
|Vo| 2 ,seeFig.[S] 



FIG. 7: (Color online) The same as Fig. gib) but for quasi-TE 
plasmon in the silicon-graphene-silicon geometry. 



B. Quasi-TE plasmon 

As follows from the dispersion relation in Eq. (|5ip . 
in order to excite TE plasmon in an assymetric geome- 
try with different substrate and cladding dielectrics one 
has to make conductivity of graphene strong enough: 
\aW\ > y/\e s — e c \. We would like to note however, 
in contrast to the case of a conventional dielectric slab 
waveguide, here only the difference between the cladding 
and substrate dielectric constants matters. Apparently, 
for the chosen doping level of graphene, TE plasmon 
does not exist in the silicon-graphene-air configuration, 
cf. Fig. H]Ja). Instead, we consider the fully symmet- 
ric configuration with silicon in the cladding and sub- 
strate: e s = e c = 12. The corresponding dispersion is 
plotted in Fig. |5] Due to the low values of |a^|, TE 
plasmon is only weakly localized. However, one benefits 
from much smaller propagation losses per plasmon period 
T = 27r/(/3fc), compared to TM plasmons. 

As the result of weak localization, typical values of 
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the effective noniinear coefficient T for TE piasmons 
are nearly twelve orders of magnitude below those for 
TM piasmons, see Fig. [7] Remarkably, for TE plas- 
mon nonlinear contribution from the dielectrics is domi- 
nant over graphene, they become comparable only when 
the frequency f2 approaches the interband absorption 
threshould ftth = 2. 

With the account of small propagation constants /3, 
diffraction term and associated effects due to its inter- 
play with the focusing nonlinearity become important for 
TE piasmons with typical widths of the order of sev- 
eral micrometres. Taking f2 = 1.9 (corresponding to 
A ~ 6.5/im), for L y = 5/xm the diffraction length be- 
comes Ld = Ly/3k « 30/im, and the plasmon period is 
27r/(/3fc) ss 0.5/im. However, due to the low nonlinearity, 
one requires considerably high powers to observe basic 
effects such as self-focusing. For instance, to form a spa- 
tial soliton of width L y the peak power density should be 
IV'oI 2 = (TLd)- 1 « 3 • 10 10 W/m, and therefore the total 
power must be of the order of several mega- Watts. 

VI. SUMMARY 

Using asymptotic expansion of Maxwell equations and 
boundary conditions, we have derived amplitude equa- 
tion for nonlinear TM and TE surface plasmon waves 
in the dielectric-graphene-dielectric planar configuration. 



Induced surface current in graphene shows strongly non- 
linear response to the applied electromagnetic field. We 
have shown that this leads to the effective focussing Kerr 
type nonlinearity. For TE plasmon this nonlinearity is 
fully analogous to that of an infitesimally thin dielectric 
layer. However, for TM plasmon the structure of the cor- 
responding nonlinear coefficient is different and reflects 
the unique surface-only response of graphene. 

For typical doping levels of graphene of the order of 
O.leV, TM piasmons are strongly localized. This causes 
the significant enhancement of nonlinearity, we predict 
that considerable nonlinear phase shifts can be observeed 
for power densities as low as few micro- Watts per metre. 
Remarkably, graphene contribution to the overall non- 
linearity is shown to be strongly dominant over that of 
dielectrics in this case. In contrast, TE piasmons are 
only weakly localized, and the major part of the overall 
nonlinearity is due to dielectric substrate and cladding. 
Typical values of the nonlinear coefficient for TE pias- 
mons are found to be about 12 orders of magnitude below 
those for TM piasmons. 
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